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ABSTRACT 

Magnetohydronamic turbulence is believed to play a crucial role in heating the laboratorial, 
space, and astrophysical plasmas. However, the precise connection between the turbulent fluc¬ 
tuations and the particle kinetics has not yet been established. Here we present clear evidence 
of plasma turbulence heating based on diagnosed wave features and proton velocity distribu¬ 
tions from solar wind measurements by the Wind spacecraft. For the first time, we can report 
the simultaneous observation of counter-propagating magnetohydrodynamic waves in the so¬ 
lar wind turbulence. Different from the traditional paradigm with counter-propagating Alfven 
waves, anti-sunward Alfven waves (AWs) are encountered by sunward slow magnetosonic waves 
(SMWs) in this new type of solar wind compressible turbulence. The counter-propagating AWs 
and SWs correspond respectively to the dominant and sub-dominant populations of the imbal¬ 
anced Elsasser variables. Nonlinear interactions between the AWs and SMWs are inferred from 
the non-orthogonality between the possible oscillation direction of one wave and the possible 
propagation direction of the other. The associated protons are revealed to exhibit bi-directional 
asymmetric beams in their velocity distributions: sunward beams appearing in short and narrow 
patterns and anti-sunward broad extended tails. It is suggested that multiple types of wave- 
particle interactions, i.e., cyclotron and Landau resonances with AWs and SMWs at kinetic 
scales, are taking place to jointly heat the protons perpendicularly and parallel. 

Subject headings: solar wind — waves — turbulence 


1. Introduction 


Turbulence is a common phenom enon in vari- 


ous environments, e.g., laboratorial (IHowes et al. 

2012 

), space (Tu & Marschl 1995: iGoldstein et al. 

1995 

: Bruno & Carbonel 2013ll. and astrophysical 
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plasmas (IZhuravleva et al. 120141) . In the helio¬ 


sphere the solar wind is deemed to be an ideal nat¬ 
ural laboratory to investigate plasma turbulence, 
since particles and fields can be measured in-situ 
and the solar wind is approximately boundless 
for the turbulent eddies/fluctuations therein. The 
magnetic fluctuations are often correlated with 
the velocity fluctuations. Together with the weak 
compressibility of the magnetic field strength, this 
indicates the domi nance of Alfvenic fluct uations 
in the turbulence ( Belcher fc DavisI Il97lh . The 
so-called Elsasser variables, which are defined as 
a sum or difference of fluctuations for velocity 
and magnetic field in Alfven unit (Z± = v ± b), 
have often been employed in theoretical calcu¬ 
lations to study the counter-propagating waves 
in Alfvenic turbulence. From the space observa¬ 
tion point of view, the nature of the minor (sub- 
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dominant) Elsasser variable remains unknown, al¬ 
though it is assumed to represent sunward/inward 
Alfven waves in the scenario of Alfvenic turbu¬ 
lence. However, this notion has not yet been con¬ 
firmed. Nevertheless, sunward propagating Alfven 
waves with the absence of anti-sunward Alfven 
waves are identified in up stream regions of the bow 
shock (jWang et al.ll2015) and in some s treams of 
the pristine solar wind ( He et al. 2015b[l . 

The minor Elsasser field may also be associated 
with c onvected magnetic structures (|Tu fc Marsch 


199,1c iBruno et al.l 120011 1 , which would also con¬ 


tribute to the obser ved excess of mag netic energy 
over kinetic energy (jChen et al.11201,11 1. Therefore, 


it is still unclear whether Zminor represents in¬ 
ward waves and/or involves mostly nonlinear in¬ 
teraction with the outward Z 


major* 


The compress¬ 
ibility of solar wind turbulence usually manifests 
itself i n pressure-balan ced structures of various 
types (|Yao et al.l 1201111 . which may come in the 
form of spaghetti-like fl ux tube stemming from 
the Sun ( Borovsky 2008h . quas i-perpendicular ki¬ 
netic slow magnetos onic waves ( Howes et al ]l2012t 
ZhBng_et_a]J_ 201^ ), or mirror-mode structures 


( Zhang et al. 200£; Yao et aDl2013ll . The nature 
of Zininor and its relation to compressibility will be 
investigated here. 

Turbulence is thought to be dissipated through 
the damping of small-scale kinetic waves or 
via currents in coherent structures, and can 
thereby he at the plasma through energy conver¬ 


sion le.g.. IVoitenko fc GoossensI l2004 iMarsch 


I 2 OO 6 I : iMarkovskii fc Vasauez 201lll . Eor exam¬ 
ple, both linear and nonlinear resonances be¬ 
tween waves and particles c an break the adia- 
batic behavior of particles (|Marsch et al.l Il982 : 


Hollweg fc IsenbergI 1200211 . Resonance diffusion 
plateaus of various types, e.g., Landau reso¬ 
nance and left/right cyclotron resonance, seem 
to be consistent with the measured contours of 
proton v elocity distributions i n fast so l ar win d 


streams ( Marsch fc Tu 2001 : iHe et al. 2015a l. 


where turbulence is dominated by Alfven waves 
from the Sun. Non-resonant (stochastic) heat¬ 
ing might be important in regions with ex¬ 
tremely lo w plasma beta, e.g., clo se to the so¬ 


lar corona ( Chandran et al. 2010ll . Dissipation 


of strong current density within coherent struc¬ 
tures m ay also be an important source of local 
heating ( Parashar et al. I l2n09t lOsman et al.ll20iTl 


Karimabadi et al. 20 lit TenBarge k, Howes 2013 : 


Wang et al. 20131: Zhang et al.l 2015b ). But how 
particles are heated in compressible MHD tur¬ 
bulence involving counter-propagating waves still 
remains to be understood. This question will also 
be studied here in detail by use of Wind in-situ 
particle and field measurements. 

2. Observations and data analysis 

The solar wind observations presented here are 
from the Wind spacecraft, and are based on the 
plasma and magnetic field data provid ed respec¬ 


tively by the 3DP and MFl instruments (|Lin et al 


_ eiy 

I 995 I : iLepping et al. Ill995l) . The time interval for 


our study is [15:00, 18:00] on 2012/12/31. Time 
profiles of parameters, i.e., proton number den¬ 
sity (A^p), bulk velocity components (V||, V±i, and 
Yl 2 ), magnetic field components (i?||, B±i, and 
B± 2 ), are illustrated in the left panels of Figure 1. 
The parallel and the first and second perpendicu¬ 
lar directions are defined as the directions of bp, 
bo X X and bp x (bg x x), where bg represents 
the mean magnetic field direction pointing anti- 
sunward in this case, and x is the GSE x-direction. 
It is surprising to find that: Vj| is positively cor¬ 
related with il||, indicating sunward SMWs prop¬ 
agating obliquely to bg; Vj_i and Vj _2 are anti¬ 
correlated respectively with B±i and indicat¬ 
ing anti-sunward propagating AWs. The finding of 
sunward propagating SMWs is also confirmed by 
the anti-correlation between A^p and Vj|. 

The right panels of Figure 1 illustrate the cross¬ 
correlation spectra calculated by means of the for¬ 


mula, CCaf, = Re(WaW/')/( 




Wb 


), where Wa 


and Wb are the wavelet coefficient spectra of a and 
b. Spectra of negative CC{Np, |B|), CC(A^p, V||), 
and positive CC(V||, H||) are clear evidence of sun¬ 
ward oblique SMWs. For oblique SMWs, the den¬ 
sity fluctuation is driven by longitudinal compres¬ 
sion/expansion, and the magnetic field pressure 
is almost in balance with the thermal pressure. 
The anti-sunward AWs are characterized by spec¬ 
tra having negative CC(Vj_, Hj,). The evidence 
of counter-propagating waves is solid, and not af¬ 
fected by the definition of bg, whether it is a global 
mean field averaged over the whole time inter¬ 
val (Figure Ib-le) or a loc al mean magnet i c field 
direction as introduced bv iHorburv et al. (20^ 
(Figure Ig-lj). 
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The Elsasser variables Z± are calculated by 
Z± = V ± Va = V ± with pq being 

the proton density averaged over the interval of 
[15:00, 18:00] UT. The trace power spectral densi¬ 
ties (PSD) of Z+, Z_, V, and Va are illustrated 
in Figure 2a and 2b. The related spectral indices 
are found by power-law fitting to be —1.55 ± 0.04, 
-1.54 ± 0.03, -1.51 ± 0.03, -1.52 ± 0.03. 


The characteristic directions of Z± fluctuations 
are estimated by apply ing the Singu l ar-Va lue- 
Decomposition method (ISantolik et al.l 120031) to 
the spectral matrix of Z±. The directions with 
maximum singular values are regarded as the 
oscillation directions. The direction with mini¬ 
mum singular value is deemed a good proxy for 
the propagation direction, if the Elsasser fluctu¬ 
ations are ideally incompressible. The spectra 
of the angles between the derived “oscillation” 
or “propagation” direction and the local mean 
magnetic field directions (6*(Z-|_,bo), 6)(Z_,bo), 
0 (k;)_,bo), 0(kk,bo)) as a function of time and 
period are illustrated in Figure 2c, 2d, 2e, 2f. The 
angle 0(Z+,bo) in Figure 2c is basically smaller 
than 30°, indicating longitudinal oscillation of Z+. 
Large angles 0(Z_,bo) in Figure 2d suggest a 
transverse oscillation of Z_. 


The oscillation difference between Z^. and Z_ 
demonstrates the intrinsically different natures of 
Z+ and Z_, with the former and latter pertain¬ 
ing to the sunward SMWs and anti-sunward AWs, 
respectively. Figure 2e shows that the majority 
of 0(k;)_,bo) is larger than 70°, while 0(k_._bo) 
in Figure 2f displays a scattered probability dis¬ 
tribution. To be cautious, 0(k^,bo) may not be 
the real propagation angle of SMWs due to the 
compressibility of Z+. The derived kl for Alfven 
waves may be either the real propagation direc¬ 
tion k_ or bg. The patches of 6*(k^, Z_) ^ 90° in 
Figure 2g may infer the nonlinear action of AWs 
(Z_) on SMWs (Z+) (Z_ • V+Z+ ^ 0). The AWs 
may be shorn by the SMWs when 0(kk, Z+) ^ 90° 
(Figure 2h), which may lead to a transverse cas¬ 
cade of AWs. 


The proton 2D velocity distribution functions 
(VDFs) in the top panels of Figure 3 illustrate the 
bi-directional asymmetric beams along the bg di¬ 
rection. Such a twofold beam pattern has, to our 
knowledge, not been reported before. As com¬ 
pared to the anti-sunward beam, which appears 
prolonged in Vji and broad in V±, the sunward 


beam looks relatively shorter in and narrower 
in Vj_. We interpret this new observational pattern 
as probable evidence of wave-particle interaction 
between protons and counter-propagating waves. 


The anti-sunward beam may be the result of 
Landau and right-cyclotron resonance with the 
outward oblique AWs at kinetic scales (general¬ 
ized kinetic Alfven waves), the wave signatures 


(Bale et al. 

20051 

Sahraoui et al.l 

2010 ; 

He et al. 

2011 

20121 : 

Podesta & GarvI 2011 

; Roberts et al. 


2013c iKlein et al.l 120141 ). The sunward beam 


may be formed via Landau resonance with ki¬ 
netic SMWs, the wave features of which have 


been studied in linear theory (|Zhao et al.l 12014 ; 
Narita fc Marsch 20151) . In other time intervals, 
without clear signatures of counter-propagating 
MHD waves, the bi-directional asymmetric beams 
in the VDFs are absent. The proton ID reduced 
VDFs in the bottom panels of Figure 3 show a 
plateau and an extended tail on the sunward 
and anti-sunward sides respectively. It is excit¬ 
ing to find that the plateau positions coincide 
well with the field-projected phase speed of quasi¬ 
perpendicular SMWs (vertical dotted line in pan¬ 
els), which indicates a Landau resonance between 
protons and sunward SMWs. The dispersion re¬ 
lation for quasi-perpendicular SMW yields ap¬ 
proximately - Vl/{l + Vl/Cl), 

with Va and Cs being the Alfven speed and 
ion sound speed respectively. The propagation 
of quasi-perpendicular SMWs can also be de¬ 
scribed by the equat ions derived for the gener - 
alized Elsasser fields ( Marsch fc Mangenevlll987 ). 
But they have never before been analyzed obser- 
vationally. SMWs in the extreme limit of fc|| <C k± 
are also named as pseudo-Alfvenic waves, which 
is passive scattered/mixed by Alfv en waves with 
k\\ <C k± I Schekochihin et al. 20091 ). Alfven waves 
with considerable A:|| may also be distorted by the 
quasi-perpendicular SMWs leading to larger k± 
due to the nonlinear term ^h^mw ’ VZj_,aw in 
the convective derivative. 


3. Summary and discussions 

On the basis of the observational findings 
presented here, a new scenario of compressible 
plasma turbulence is proposed and illustrated 
in Figure 4. Dominant AWs are propagating 
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anti-sunward and quasi-perpendicular SMWs are 
propagating upstream against the AWs. The 
nonlinear interaction during collisions between 
these counter-propagating waves may lead for 
both waves to a cascade of wave energy towards 
plasma kinetic scales. Meanwhile, in response 
to the counter-propagating waves, bi-directional 
asymmetric beams emerge in the proton VDFs 
through resonant wave-particle interactions. In 
this way, the turbulence energy cascading and 
dissipation as well as the solar wind ion heating 
are included in this common scenario. This new 
type of turbulence, which is unusual in the solar 
wind fast streams, is different from the traditional 
paradigm, where counter-propagating waves are 
both Alfvenic but not yet confirmed observation- 
ally. 


Some new challenges are raised along with these 
findings. (1) How are the sunward SMWs gener¬ 
ated in the solar wind turbulence: Is it by com¬ 
pression be tween streams, the par ametric decay 
from AWs ( Verscharen et all 2012 1. or by other 
processes? (2) The observed compressible MHD 
turbulence and the concurrent collisionless damp¬ 
ing of kinetic waves, together call for a more 
comprehensive study, involving numerical mod¬ 
eling and observational statistical analysis. (3) 
What is the radial evolution in the solar wind of 
the compressible turbulence, involving AWs and 
SMWs, which have also been identified in the 
solar wind source region? (4) Besides the lo¬ 
cal wave-particle interactions, effects of large-scale 
evolution (e.g., anisotropic cooling during expan¬ 
sion, gravity, charge-separation electric field, mir¬ 
ror force, and differential flows between species) 
need to be t aken into account for the proton ki¬ 
netic s (e.g., Cranmei 2014 : Isenberg fc Vasau^ 


2015 : Hellinger fc Travnfcek 2015ll . These are is- 
sues that may well be addressed by the upcom¬ 
ing missions Solar Orbiter and Solar Probe Plus, 
which will be ’’sailing against the solar wind” into 
the inner heliosphere and outer corona. 
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Fig. 1.— Evidence of sunward quasi-perpendicular slow magnetosonic waves (SMWs) coexisting with anti- 
sunward Alfven waves (AWs); the evidence is based on the correlation analysis of various parameter pairs. 
Time profiles of A^p (black) and |B| (red) in (a), A^p (black) and V|| (red) in (b), V|| (black) and (red) 
in (c), V±^i (black) and B±^i (red) in (d), V ±^2 (black) and (red) in (e). Correlation spectra between 
different pairs of parameters (A^p, |B|) in panel (f), (A^p, V||) in (g), (Vn, i3||) in (h), (V(l,i, B±,i) in (i), 
and (Vj_. 2 , B±, 2 ) in (j). The indices T, 1 and T, 2 refer to the directions of bg x x and bg x (bg x x). The 
direction x represents the GSE-x direction. The direction bg in the left and right panels refers to the global 
and local mean magnetic field direction, respectively. 
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Fig. 2.— The Elsasser variables reveal the fluctuations to be longitudinal oscillations of sunward quasi¬ 
perpendicular SMWs for Z+ and transverse oscillations of anti-sunward AWs for Z^. Trace power spectra 
of Z+ & Z_ in (a), and V & Va in (b) during the time interval, with their fitted spectral indices being 
— 1.55 ±0.04, —1.54 ±0.03, —1.51 ±0.03 and —1.52 ±0.03, respectively, (c) Spectrum of angle between Z+ 
and bo, where Z+ direction is the characteristic direction with maximum singular values as obtained from 
Singular Value Decomposition method, and bg direction is the local mean magnetic field direction. Small 
angles coded in blue indicate the longitudinal oscillation of Z+ along bg. (d) Spectrum of 6*(Z_,bo) with 
large angles in red showing oscillation of Z_ transverse to bg. (e) Spectrum of 0(k^, bg), where k^J, direction 
that is basically representing the eigenvector with minimum variance of Z_|_ may have some deviation from 
the real propagation direction due to the compressibility, (f) Spectrum of 6*(kl,bg) with kl being the 
characteristic direction of Z_ minimum variance. The possible nonlinear impact of Z_ on Z+ may be 
inferred from the patches with 6*(k^,Z_) < 90° in (g). Patches with small angles in (h) imply the possible 
role of SMWs in cascading the energy of outward AV^. 
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Fig. 3.— Bi-directional asymmetric proton beams associated with the counter-propagating waves, indicating 
the solar wind heating by turbulence via Landau and cyclotron resonances, (al) Cut view of proton velocity 
distributions at 16:11:23 with the GSE-x direction (R* defined as -x in GSE) and magnetic field direction 
(Bo ) in the plane, showing sunward and anti-sunward beams on either side of the core are aligned with the 
magnetic field direction (Bq, white line). Different from the anti-sunward beam, which extends further from 
the core and has a broad profile across Bq, the sunward beam drifts at a smaller speed and appears to have 
a narrower transverse velocity profile, (bl) Reduced VDF as integrated from the VDF in the plane of (Bq, 
R* X Bq) along the direction of R* x Bq. Left dotted line indicate the position of —Cg (the phase speed of 
quasi-perpendicular slow mode wave as projected along Bq), and right dashed line is located at the position 
of Va (Alfven speed). Velocity distributions at another three times are illustrated in (a2, b2), (a3, b3), and 
(a4, b4). 
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Fig. 4.— The sketch illustrates solar wind turbulence that is characterized by collisions between outward 
AWs and inward SMWs. These collisions may lead to wave coupling, energy cascading, dissipation and 
heating. Transverse oscillations of Z_ for outward AWs are described with three wave forms along three 
background magnetic field lines. Longitudinal oscillation of for inward SMWs are represented by two 
wave forms aligned perpendicular to Bq. The trajectory of the WIND spacecraft in the solar wind frame is 
denoted by the green dashed line. 
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